SELF-SIMILAR SOLUTIONS OF THE NON-STRICTLY HYPERBOLIC 

WHITHAM EQUATIONS 

V. U. PIERCE AND FEI-RAN TIAN 

Abstract. We study the Whitham equations for the fifth order KdV equation. The equations are 

o , 

\ neither strictly hyperbolic nor genuinely nonlinear. We are interested in the solution of the Whitham 

equations when the initial values are given by a step function. We classify the step like initial data 
C***'" ■ into eight different types. We construct self-similar solutions for each type. 

(N 

£ ■ 1. Introduction 

It is known that the solution of the KdV equation 

(1.1) u t + 6uu x + e 2 u xxx = 

O 

has a weak limit as e — > while the initial values 

O 

• S ■ u(x,0;e) = u Q (x) 
> 

•rH . are fixed. This weak limit satisfies the Burgers equation 

X 

(1.2) ut + (3u 2 ) x = 

until its solution develops shocks. Immediately after, the weak limit is governed by the Whitham 
equations H C2 031 

(1.3) u it + Ai(ui,u 2 ,U3)wfa = , * = 1, 2, 3, 

where the Aj's are given by formulae ()2.12j) . After the breaking of the solution of (jl.3|) . the weak 

limit is described by systems of at least five hyperbolic equations similar to (jl.3j) . 

l 
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The KdV equation is just the first member of an infinite sequence of equations, the second 
of which is the so-called fifth order KdV equation 

(1.4) u t + 30m u x + 20e u x u xx + lOe uu xxx + e u xxxxx = . 

The solution of the fifth order KdV equation (jl.4|) also has a weak limit as e — ► 0. As in the KdV 
case, this weak limit satisfies the Burgers type equation 

(1.5) m 4 + (10u 3 ) x = 

until the solution of (jl.5j) forms a shock. Later, the limit is governed by equations similar to (jl.3j) . 
namely, 

(1.6) u it + Hi(ui,u 2 ,u 3 )u ix = , i = l,2,3, 

where /Vs are given in ([2. 18)1 . They will be also be called the Whitham equations. 

In this paper, we are interested in the solution of the Whitham equation (jl.6|) for the fifth order 
KdV fjl.4|) with a step like initial function 

{a x < 
a ^ b . 
b x>0 

For such an initial function with a>0, b<a or a<0,b> a, the solution of the Burgers type 
equation (J1.5)) has already developed a shock at the initial time, t = 0. Hence, immediately after 
t = 0, the Whitham equations (|1.6J) kick in. Solutions of (jl.6j) occupy some domains of the space-time 
while solutions of ()1.5|) occupy other domains. These solutions are matched on the boundaries of the 
domains. 

The solution of the Burgers equation (jl.2j) with initial function (|1.7J) is simple: it is either a 
rarefaction wave or a single shock wave. 

The Burgers type equation (jl.5|) is more complicated, as its flux function changes convexity at 
u = 0. Its solution with step like initial data (jl.7J) can be a rarefaction wave, a single shock wave 
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or a combination of both jE] . As a consequence, the solutions of the Whitham equations p. 6)1 will 
be seen to be more complex than those of p. 3)1 in the KdV case. Indeed, there are eight types of 
different solutions in the former case while there is only one type of solution in the latter case. 
The KdV case with the step like initial data (jl.7[) was first studied by Gurevich and Pitaevskii 
The Burgers solution of (jl.2)l develops a shock only for a > b. Moreover, the corresponding 
initial function is equivalent to the case a = 1, b = 0. In this case, Gurevich and Pitaevskii found 
a self-similar solution of the Whitham equations p.3jl . Namely, the space-time is divided into three 
parts 

(1) |<-6, (2) "6<y<4, (3) |>4. 
The solution of (jl.2J) occupies the first and third parts, 

(1.8) u(x, t) = 1 when ~ < — 6 , u(x,t) = when f > 4 . 

The Whitham solution of p. 3)1 lives in the second part, 

x 

(1.9) ui(x,t) = l, - = A 2 (l, «2,0) , u 3 (x,t) = 0, 
when — 6 < x/t < 4. 

Whether the second equation of (jl.9|) can be inverted to give U2 as a function of the self-similarity 
variable x/t hinges on whether 

0M % 1*2,0) ^0. 



Indeed, Levermore has proved the genuine nonlinearity of the Whitham equations (|1.3|) . i.e., 
(1.10) ^K^.ua) >0, i = 1,2,3, 

OUi 

for u\> U2> U3. 

For the fifth order KdV (jl.4J) . equations p. 6)1 . in general, are not genuinely nonlinear, i.e., a 

property like p. 10)1 is not available. Hence, solutions like p. 8)1 and p. 9)1 need to be modified. 
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Our construction of solutions of the Whitham equation ()1.6|) makes use of the non-strict hyperbol- 
icity of the equations. For KdV, it is known that the Whitham equations (jl.3j) are strictly hyperbolic, 
namely: 

X 1 (u 1 ,u 2 ,u 3 ) > \ 2 {u 1 ,u 2 ,u 3 ) > X 3 (u 1 ,u 2 ,u 3 ) 

for ui > u 2 > u 3 . For the fifth order KdV ()1.4j) . different eigenspeeds of (jl.6j) . fii(ui,u 2 ,u 3 )'s, may 
coalesce in the region u x > u 2 > u 3 . 

For the fifth KdV with step-like initial function (jl.7|) where a = 1 and 6 = 0, the space time is 
divided into four regions (see Figure 1.) 

(1) | < -15 , (2) - 15 < - < a , (3) a < | < 16 , (4) | > 16 , 

where a is determined by ()3.15|) . In the first and fourth regions, the solution of (|1.5|) governs the 
evolution: 

u(x,t) = 1 where x/t < —15 and u(x,t) = where x/t > 16 . 

The Whitham solution of (jl.6)l lives in the second and third regions; namely: 

x x 
(1.11) ui(x,t) = l, - = fi 2 (l,U 2 ,U 3 ) , - = fl 3 (l,U 2 ,U 3 ) , 

when —15 < x/t < a, and 

x 

ui(x,t) = l, - = /i 2 (l, «2,0) , u 3 (x,t)=0, 

when a < x/t < 16. 
Equations ()1.11|) yield 

j4 2 (l, u 2 ,u 3 ) = fi 3 (l, u 2 ,u 3 ) 

on a curve in the region < u 3 < u 2 < 1. This implies the non-strict hyperbolicity of the Whitham 
equations (J1.6)) for the fifth order KdV. 

The organization of the paper is as follows. In Section 2, we will study the eigenspeeds, /ij(wi, u 2 , u 3 ) 's, 
of the Whitham equations (jl.6j) . In Section 3, we will construct the self-similar solution of the 
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Figure 1. Self-Similar solution of the Whitham equations for a = 1 and b = of type II. 

Whitham equations for the initial function ()1.7j) with a = 1, b = 0. In Section 4, we will use the 
self-similar solution of Section 3 to construct the minimizer of a variational problem for the zero 
dispersion limit of the fifth order KdV. In Section 5, we will consider all the other possible step like 
initial data ()1.7|) . We find that there are eight different types of initial data. We construct self-similar 
solutions for each type. 

2. The Whitham Equations 

In this section we define the eigenspeeds of the Whitham equations for both the KdV (jl.lj) and 
fifth order KdV (jl.4|) . We first introduce the polynomials of £ for n — 0, 1, 2, . . . 13*1 HUj: 



(2.1] 



P n (£, ui, u 2 , u 3 ) = C +1 + a n ,iC H h a>n,n+\ 



where the coefficients, a nj i, a nj2 , . . . , a n) „ + i are uniquely determined by the two conditions 

P n (£,ux,u 2 ,u 3 ) 



(2.2) 
and 
(2.3) 



C' 1/2 + C(r 3/2 ) for large |f| 



"2 



P„(£, m, u 2 , u 3 ) 



"3 



Here the sign of the square root is given by vT = 1- 
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In particular, 

(2.4) P (^,u 1 ,u 2 ,u 3 ) = £ + a ,i , Pi(^,u 1 ,u 2 ,u 3 ) = £ 2 - -(u x + u 2 + u 3 )£,+ai t2 

where 

/ Ms) 

K(s) 

«1,2 = \{u\U 2 + M1M3 + M 2 n 3 ) + -{ui +u 2 + u 3 )a 01 . 
6 6 



Here 



u 2 - u 3 
s = 



Ui - u 3 

and K(s) and E(s) are complete elliptic integrals of the first and second kind. 
K(s) and E(s) have some well-known properties |Hl I0J - As — 1 < s < 1, we have 

7r r s 9 o , 1 • 3 • • • (2ra — 1) N r, „ 

2.5 K(s) = - 1 + - + — s 2 + ■■■ + ( i Lf s n + ... 

v 7 v ; 2 1 4 64 v 2 • 4 • ■ ■ 2n ' J ' 

(2.6) = ^[l-£-A S 2 1_ l-3---(2n-l) ^ 

v ; v ; 2 L 4 64 2n - l v 2 • 4 • • • 2n ' 

while, as 1 — s <C 1, we have 

1 16 

(2.7) K(s) » -log 



1 16 

(2.8) £(s) « l + -(l- S )[log- 1] . 



2 1 - s 

4 V /L 1 - s 



Furthermore, 
(2.9) 
(2.10) 



dK(s) E(s) - (1 - s)K(s) 



ds 2s(l - s) 

d£(s) £(s) - K(s) 



(is 2s 
It immediately follows from ()2.5|) and ()2.6j) that 

(2,1) 1 <|W<^i /or0 < S <l 
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The eigenspeeds of the Whitham equations (jl.3j) are defined in terms of P and Pi of (|2.4J1 . 

X i (u 1 ,u 2 ,u 3 ) = 12— -, z = l,2,3, 

P (u i ,u 1 ,u 2 ,u 3 ) 

which give 

K(s) 



Xi(u 1 ,u 2 ,u 3 ) = 2(m + u 2 + u 3 ) +i(ui -u 2 ) 
(2.12) A 2 (wi, u 2 , u 3 ) = 2(u l + u 2 + u 3 ) + A(u 2 - u x ) 



E(s) ' 

sK(s) 



E(s) - (1 - s)lf(s) ' 
K(s) 



X 3 (u u u 2 ,u 3 ) = 2(u 1 + u 2 + u 3 )+A(u 2 -u 3 )- 

E(s) — K[s) 

Using (|2.1ip . we obtain 

(2.13) Ai -2( Ul + u 2 + u 3 ) > 0, 

(2.14) A 2 -2( Ul + u 2 + u 3 ) < 0, 

(2.15) A 3 -2{ Ul + u 2 + u 3 ) < 0, 

for u\ > u 2 > u 3 . In view of (j2.5M2.8|) . we find that Ai, X 2 and A3 have behavior: 
(1) At u 2 = u 3 : 



Xi(u-i_,u 2 ,u 3 ) = 6«i , 

A 2 (ui, u 2 , u 3 ) = X 3 (ux,u 2 , u 3 ) = 12u 3 - 6ui . 



(2.16) 

(2) At ui = u 2 : 

Ai(«i, u 2 ,u 3 ) = X 2 (u!,u 2 ,u 3 ) = 4ux + 2u 3 , 

(2.17) 

A 3 (iti, 1*2,1*3) = 6% . 
The eigenspeeds of the Whitham equations (jl.6|) are 

/ q1q n / \ Qn P2(u i ,U 1 ,U 2 ,U 3 ) . 

(2.18) //i(iii, m 2 , u 3 ) = 80— , 1 = 1, 2, 3 

ii)(«»,«l,«2,«3j 

They can be expressed in terms of Ai, A 2 and A 3 of the KdV. 
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Lemma 2.1. [TU] The eigenspeeds, fii's, satisfy: 



1 dq 
(2.19) u 2, W3) = 7: [A, - 2(ni + u 2 + w 3 ) 77— «2, M3) + q(ui, «2, W3) , 

where q(ui, 112,113) is the solution of the boundary value problem of the Euler-Poisson-Darboux 
equations: 



g(it, w, it) = 30m 



2 



(2.21) g = _^_ frl( _ tt] , MJ . 

The solution of (|2.2(Jj) is a symmetric quadratic function of ui, u 2 and u 3 

(2.22) 5(^1, u 2 , W3) = 6(«i + «2 + w 3 ) + 4(ui« 2 + W1U3 + W2M3) 
For KdV, Ai's satisfy 

<9A 3 3 A 2 - A 3 <9A 2 

<9« 3 2 W 2 — W 3 <9w 2 

for W3 < u 2 < U\. Similar results also hold for the fifth order KdV. 
Lemma 2.2. 

(2.23) p. > Vp->0, 

ou 2 2 u 2 — u 3 ou 2 

(2.24) p < 3 -^l ifP->0, 
for u 3 < u 2 < u\ . 
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Proof. We use (|2.19|) to calculate 

<9/i 3 1<9A 3 dq 1 <9 2 g 
^ = 2^^ + 2 [A3 - 2(Ml+M2 + M3)] ^| 

(2-25) < |L + J [As - 2( Ml + « 2 + W3 )]f| , 

4« 2 - w 3 aw 3 2 c/Mg 



and 



1 <9(/ 1 / 
A*2 - A*3 = ~ (A2 - A 3 ) h -[A 3 - 2(ui + w 2 + M3) 1 



2 <9w 3 2 " \^ M 2 

1 dq 1 <9 2 <7 

77 (A 2 - A 3 ) -7^ + o [A3 - 2(ui + w 2 + w 3 )]2(w 2 - u 3 ) 



2 <9w 3 2 du 2 du 3 

2, , / 3 A 2 — A 3 dq 3 ri , d 2 q 
(2.26) = -(« 2 -« 3 ) 7— -7 i + o[ A 3- 2 («i + «2 + ^]- 



3 \4m 2 - u 3 du 3 2 du 2 du 3 

where we have used equation (|2.2U|) 

dq dq d 2 q 
2(u 2 - u 3 ) 



du 2 du 3 du 2 du 3 
It follows from formula ()2.22j) for q that 

d 2 q d 2 q 

du 2 du 3 du\ 

which, along with with (jZSHD and (f2~2TTj) . proves ()2.23j) . 
Inequality (j2.24j) can be proved in the same way. 



The following calculations are useful in the subsequent sections. 

Using formula (j2.19J) for \i 2 and ^3 and formulae (j2.12|) for A 2 and A3, we obtain 

2(w 2 — u 3 )K 

(2.27) /ia(wi, «2, U3) - A*3(wi, «2, "3) = tj^ _ - (1 - 9) A - ] M2 ' 



where 



M (ui , U2 , B3 , = ^ + ( 1 -^-( 1 -< + ^. 
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We then use (EH), dUEj) and (jHHD to calculate 

/ 2 2g N aM(ni,M 2 ,w 3 ) = 10K - 3n 2 - m 3 ) _ 

<9u 2 «i - «3 

We next consider 

(2.29) F(u 1 ,u 2 ,u 3 ) := • 

M 2 - W 3 

Using formula (|2.19|) for // 2 and /z 3 and formulae (|2.12|) for A 2 and A 3 , we obtain 
F = 2 (1 ~ S)K 9q I 2 * 9g 



£ — (1 — s)i^ <9w 2 K-Edu 3 

s(l-s)K _ d 2 g K (l-s)K dq 

E-(l-s)K {Ul U3) du 2 du 3 + [ K-E E-(l-s)K*du 3 ' 

where we have used equations (|2.2(J|) in the last equality. Finally, we use the expansions (j2.5H2.6p for 

K and E to obtain 

7 <9 2 <? 3 2 <9g . 



F( Ul ,« 2)W3 ) = -4[(2 --« + •..)(«!- lis) + (- 7 + 0(^)) 



4 du 2 du 3 4 <9ii 3 

(2.30) -16[(2 - ~s + •••)(«! - + (~ + 0(s 2 ))( Wl + « 2 + 3« 3 )] , 

where we have used formula ()2.22|) for q in the last equality. 

3. A Self-similar Solution 

In this section, we construct the self-similar solution of the Whitham equations (|1.6|) for the initial 
function (jl.7|) with a = 1 and 6 = 0. We will study all the other step like initial data in Section 5. 

Theorem 3.1. (see Figure 1.) For the step-like initial data Uq(x) of Ji.7| ) with a = 1,6 = 0, the 
solution of the Whitham equations M.b\) is given by 

(3.1) Mi = 1 , X = /i 2 (l,M 2 ,M 3 ) t , x = fx 3 (l,u 2 ,u 3 ) t 
for — 15£ < x < at and by 

(3.2) u\ = 1 , x = /i 2 (l,w 2 ,0) t , u 3 = 
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for at < x < 16t, where a = /^(l, u * 0) and u* is the unique solution u 2 of ^t 2 (lj u 2, 0) = /^(l, u 2 , 0) 
in the interval < u 2 < 1. Outside the region — 15t < a; < 16i, £/ie solution of the Burgers type 
equation Ji.,5)) zs g'wen 5y 

(3.3) u = 1 x < -15t 
and 

(3.4) u = x > I6t . 

The boundaries x = —15t and x = 16t are called the trailing and leading edges, respectively. They 
separate the solutions of the Whitham equations and Burgers type equations. The Whitham solution 
matches the Burgers type solution in the following fashion (see Figure 1.): 



(3.5) ui = the Burgers type solution defined outside the region , 

(3.6) u 2 = u 3 , 

at the trailing edge; 

(3.7) ui = u 2 , 

(3.8) u 3 = the Burgers type solution defined outside the region , 
at the leading edge. 



The proof of Theorem 3.1 is based on a series of lemmas. 

We first show that the solutions defined by formulae (j3.1|) and (j3.2j) indeed satisfy the Whitham 
equations (TEJ) [HE]. 

Lemma 3.2. (1) The functions u\, u 2 and u 3 determined by equations \3. 1)) give a solution of 
the Whitham equations M.b]) as long as u 2 and u 3 can be solved from as functions of x 
and t. 
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(2) The functions U\, u 2 an d u 3 determined by equations h3. 2]) give a solution of the Whitham 
equations M.n\) as long as U2 can be solved from \H. jjj) as a function of x and t. 

Proof. (1) u\ obviously satisfies the first equation of (jl.fij) . To verify the second and third equations, 
we observe that 

(3.9) ^ = ^ = 

du* duo 



on the solution of (|3.1jl . To see this, we use (|2.21J) to calculate 



9^2 _ du 3 



d\ 2 

(A*2 - ^3) = . 



<9w 3 A 2 - A 2 

The second part of (J3.9j) can be shown in the same way 

We then calculate the partial derivatives of the second equation of (|3.1|) with respect to x and t. 

d/i 2 diii 

1 = tU 2x , = -— tu 2t + [ll , 

ou 2 ou 2 

which give the second equation of (jl.6|) . 

The third equation of p.fjjl can be verified in the same way. 
(2) The second part of Lemma 3.2 can easily be proved. 

□ 

We now determine the trailing edge. Eliminating x and t from the last two equations of (|3.1|) 
yields 

(3.10) /i 2 (l, u 2 , u 3 ) - ji 3 (l, u 2 , u 3 ) = . 
Since it degenerates at u 2 = w 3 , we replace (j3.1U|) by 

(3.11) F(l,u 2 ,u 3 ) := = 0. 

u 2 - u 3 

Here, the function F is also defined in ()2.29|) . 
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Therefore, at the trailing edge where w 2 = w 3 , i.e., s = 0, equation (|3.11|1 . in view of the expansion 
(I2.30J) . becomes 

2(l-w 3 )-|(l + 4« 3 ) = 0, 

which gives u 2 = u% = 1/4. 

Lemma 3.3. Equation has a unique solution satisfying 112 = 113. The solution is w 2 = w 3 = 

1/4. The rest of equations \3. 1\) at the trailing edge are u\ = 1 and x/t = ^ 2 (1, 1/4, 1/4) = —15. 

Having located the trailing edge, we now solve equations (|3.1|) in the neighborhood of the trailing 
edge. We first consider equation (|3.11|) . We use ()2.30|) to differentiate F at the trailing edge U\ = 1, 
u 2 = u 3 = 1/4 



40, 



du 2 du 3 

which show that equation (j3.11j) or equivalently (|3.1()jl can be inverted to give M3 as a decreasing 
function of U2 

(3.12) u 3 = A(u 2 ) 

in a neighborhood of m 2 = m 3 = 1/4. 

We now extend the solution ()3.12|) of equation (|3.10j) in the region 1 > m 2 > 1/4 > m 3 > as far 
as possible. We deduce from Lemma 2.2 that 

(3.13) ^>0, ^<0 

OU2 ou 3 

on the solution of ()3.10|) . Because of ()3.9|) and ()3.13|) . solution ()3.12|) of equation ()3.10|) can be 
extended as long as 1 > u 2 > 1/4 > u 3 > 0. 

There are two possibilities: (1) m 2 touches 1 before or simultaneously as m 3 reaches and (2) m 3 
touches before w 2 reaches 1. 

It follows from (j2~T7jl and (J2~mj) that 

/x 2 (l, 1, u 3 ) > /i 3 (l, l,tt 3 ) for < m 3 < 1 . 
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This shows that (1) is impossible. Hence, U3 will touch before u 2 reaches 1. When this happens, 
equation ()3.1()j) becomes 

(3.14) ^ 2 (l,u 2> 0)-// 3 (l,w a ,0) = 0. 

Lemma 3.4. Equation \3. 14\ ) has a simple zero in the region < u 2 < 1, counting multiplicities. 
Denoting the zero by u* , then /x 2 (l, u 2 , 0)— 7x3(1, u 2 , 0) is positive for u 2 > u* and negative for u 2 < u* . 

Proof. We now use (|2.27|) and ()2.28|) to prove the lemma. In equation (|2.27|) . K — E and E — (1 — s)K 
are all positive for < s < 1 in view of (j2.11j) . By ()2.28j) . 

dM(l,u 2 ,0) = for0<« 2 <l. 
du 2 

Since M(l, m 2 , 0) vanishes at u 2 = and is positive at u 2 — 1 in view of (j2.5H2.8|) . we conclude from 
the above derivative that M(l,u 2 ,0) has a simple zero in < u 2 < 1. This zero is exactly w* and 
the rest of the theorem can be proved easily. □ 

Having solved equation ()3.10|) for u 3 as a decreasing function of u 2 for 1/4 < u 2 < u*, we turn to 
equations (|3.1|) . Because of ()3.13|) . the second equation of ()3.1|) gives u 2 as a increasing function of 
x/t, for —15 < x/t < a, where 

(3.15) a = /x 2 (l,u*,0). 
Consequently, M3 is a decreasing function of x/t in the same interval. 

Lemma 3.5. The last two equations of can be inverted to give u 2 and u% as increasing and 

decreasing functions, respectively, of the self- similarity variable x/t in the interval —15 < x/t < a, 
where a = /i 2 (l,w*,0) and u* is given in Lemma 3.4- 

We now turn to equations ()3.2j) . We want to solve the second equation when x/t > a or equivalently 
when u 2 > u*. According to Lemma 3.4, /x 2 (l, u 2 , 0)— /x 3 (l, u 2 , 0) > for u* < u 2 < 1, which, together 
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with ()2.23j) . shows that 

^ 2 (1,^2,0) 
du 2 

Hence, the second equation of ()3.2|) can be solved for u 2 as an increasing function of x/t as long as 
u* < 1*2 < 1- When u 2 reaches 1, we have 

x/t = /i 2 (l, 1, 0) = 16 , 

where we have used ()2.17|) and (|2.19|) in the last equality. We have therefore proved the following 
result. 

Lemma 3.6. The second equation of can be inverted to give u 2 as an increasing function of 
x/t in the interval a < x/t < 16. 

We are ready to conclude the proof of Theorem 3.1. 
The Burgers type solutions ()3.3j) and ()3.4|) are trivial. 

According to Lemma 3.5, the last two equations of ()3.1)1 determine u 2 and 113 as functions of x/t 
in the region —15 < x/t < a. By the first part of Lemma 3.2, the resulting u±, u 2 and U3 satisfy the 
Whitham equations (jl.6|) . Furthermore, the boundary conditions (|3.5|) and (J3.6)) are satisfied at the 
trailing edge x = —15t. 

Similarly, by Lemma 3.6, the second equation of ()3.2|) determines u 2 as a function of x/t in the 
region a. < x/t < 16. It then follows from the second part of Lemma 3.2 that u\, u 2 and U3 of ()3.2j) 
satisfy the Whitham equations (jl.6j) . They also satisfy the boundary conditions ()3.7|) and (j3.8J) at 
the leading edge x = 16t. 

We have therefore completed the proof of Theorem 3.1. 

A graph of the numerical solution of the Whitham equations is given in Figure H 
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4. The Minimization Problem 

The zero dispersion limit of the solution of the fifth order KdV equation (jl.4|) with step-like initial 
function (|1.7|) . a = 1, b = 0, is also determined by a minimization problem with constraints jUIHlE] 

1 f 1 f 1 7? — a f 1 

(4.1) Minimize! / / log ipMipi^dr/dfi + / [rjx — 16r] 5 t]?p(r))dri} . 

{^>o, ^eL 1 } 2tt Jo J q rj + fi J 

In this section, we will use the self-similar solution of Section 3 to construct the minimizer. We 

first define a linear operator 

2tt Jo \V + Vj 

The variational conditions are 

(4.2) Lip = xr] — 16tr] 5 where ip > , 

(4.3) Lip < xi] — I6trf where ip — . 
The constraint for the minimization problem is 

(4.4) ip > . 
The minimizer of (14.1)1 is given explicitly: 

Theorem 4.1. The minimizer of the variational problem \4-l\) is as follows: 
(1) Forx < -15t, 

—xi] + 80t7](r] 4 — \tf — ~) 

W) = 7= 7 — • 

V 1 - T 

(%) For -15t < x < at, 



xrjPp (ri 2 , 1 ,u 2 ,U3 )+80triP 2 (t? 2 , 1 ,u 2 ,u 3 ) ^ /^j- 

V (l~»7 2 )(«2-») 2 )(«3-'7 2 ) X 



^(77) = < 



v /% < 1] < ^/U~ 2 , 

-rriyPo (rj 2 ,l,u 2 ,0)+80tr ) P 2 (rj 2 , 1 ,u 2 ,0) 



\/ (l-»? 2 )('? 2 -'"2)(r? 2 -«3) 



U 2 < Tj < 1 
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where Pq and P<i are defined in \2.1\l and and u% are determined by equations \3.1\) . 
(3) For at<x< IQt, 



Tj < y/U^ 

xP (r} 2 ,l,u 2 ,0)+80tP 2 (r] 2 ,l,u 2 ,0) 



U2 < V < 1 



where u 2 is determined by \3. fy) . 
(4) Forx > 16t, 

t/j(rj) = . 

Proof. We extend the function ip denned on [0, 1] to the entire real line by setting ip{v) = f° r V > 1 
and taking ip to be odd. In this way, the operator L is connected to the Hilbert transform H on the 
real line [3]: 



Lij)(r}) 



P 1 f +co ib(u) 

/ Hi/)(fi)dfj, where H^(rj) = -P.V. / -t-^-dfi . 

JO * J-oo V-f* 



We verify case (4) first. Clearly ip(rj) = satisfies the constraints (|4.4j) . We now check the 
variational conditions (|4.2H4.3|) . Since ij) — 0, 

i/0 = < a;r/ — 16t/7 5 , 

where the inequality follows from x > 16t and < 77 < 1. Hence, variational conditions f|4.2H4.3|) are 
satisfied. 

Next we consider case (1). We write as the real part of (71(77) for real 77, where 



01 = y— 1 x — 80tr? H ; — — . 

The function g\ is analytic in the upper half complex plane Im(rj) > and 171(77) ~ 0(1 /r) 2 ) for large 
\tj\. Hence, Hip(rj) = Im[gi(rf)\ = x — 80trj A on < 77 < 1, where H is the Hilbert transform Q . We 
then have for < 77 < 1 



L?/>(77) 



rv 

j Hr)(ii)dfjL = xt] — 16tr] 5 
Jo 
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2 -!) 2 >o, 



which shows that the variational conditions are satisfied. Since 15+80(?7 4 — |r/ 2 — |) = 80(77 
it follows from x < — 15i that if> > 0. Hence, the constraint (|4.4j) is verified. 

We now turn to case (2). By Lemma 3.5, the last two equations of (j3.1j) determine u 2 and u 3 as 
functions of the self-similarity variable x/t in the interval —15 < x/t < a. 

We write if) = Re (g 2 ) for real 77, where 



92 



-l(x - 80tr] 4 



^ll-xriPoir] 2 , 1, u 2 , u 3 ) + 80tr]P 2 (r] 2 , 1, u 2 , u 3 )] 



y/r) 2 - l)(rf - u 2 )(r] 2 - u 3 ) 

The function g 2 is analytic in Im(r)) > and g 2 {rj) ~ 0(l/r] 2 ) for large \rj\ in view of the asymptotics 
(|2.2j) for Pq and P 2 . Hence, taking the imaginary part of g 2 yields 



Hip(r}) = < x - 80trj 4 - [ - 



x - 80tr/ 4 

-xP (r ? 2 ,l,M2,0)+80tP 2 (?? 2 ,l,U2,0)]r? 



V (l-»? 2 )(«2-r) 2 )0) 2 -«3) 

x - 80t/7 4 



< 77 < ^u~i 

Ui < 77 < y/U2~ ■ 
U2 < Tj < 1 



We then have 



xt] — IQtrf 



0<T]< ^ 



Lip{rj) = < 



XT] — 16t77 5 — J . 



"3 ^/(l-^)^-^)^-^) 
XT] — 16t77 5 



■U2 < ^ < 1 



where we have used 



(4.5) 



/f 2 



[-xP + 80tP 2 ]/j 



-dfi = 



a/(1 - H 2 )(u 2 - /i 2 )(/i 2 - m 3 ) 

which is a consequence of (|2.3|) for P an d ^2- 

We study the zeros of — xPq + 80tP2. It has two zeros at rj = ^fu^ and 77 = ■Ju 3 . This follows 
from (|2.18|) and (|3.1j) . It also has a zero between ^fu^ and ^Jui because of (j4.5j) . Since it is a cubic 
polynomial of ?7 2 , — xPo + 80tP2 has no more than three zeros on the positive 77 axis and furthermore 
these three positive zeros are simple. 
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Since the leading term in — xPq + 80tP 2 is 80t?7 6 , the polynomial is positive for 77 > ^fu^ and 
negative for < 77 < y/ui. This proves if) > 0; so (|4.4|) is verified. Since — xP + 80tP 2 changes sign 
at each simple zero, it follows from ()4.5j) that 

[-xP + 80tP 2 ]/i 

=q/i > 
»3 v(! - /" 2 )(«2 - M 2 )(/J 2 - u 3 ) 



for y/u3 < i] < \fu2. This verifies the variational conditions (|4.2j) and f|4.3|) . 

We finally consider case (3). By Lemma 3.6, the second equation of ()3.2|) determines u 2 as an 
increasing function of x/t in the interval a < x/t < 16. 

We write ip = Re(g%) for real 77, where 

y/(r] 2 - l)(r] 2 - u 2 ) 

The function g% is analytic in Im(r]) > and (73(77) ~ 0{l/rf) for large |?7| in view of the asymptotics 
()2.2j) for Po and P 2 . Hence, taking the imaginary part of g% yields 

L _ 8(^4 _ -xPp(,M , 2 ,0) + 8Q f P 2 (,M,n 2 ,0) Q<r/< 

I x - 80t?7 4 v /u 2 " < 77 < 1 

We then have 

( xrj - mtrj 5 - r xP ^ mtp \ dv o< v <^r 2 

I xr] — 16trj 5 y/u^ < 77 < 1 

where we have used 

f46) -xP (/x 2 , 1 tx 2 , 0) + 80tP 2 (/i 2 , 1, tx 2 , 0) d =q 

which is a consequence of (12.3)1 for P and P 2 . 

The function — xP (t7 2 , 1, u 2 , 0) + 80/;P 2 (77 2 , 1, u 2 , 0) has two zeros on the positive 77-axis. One is at 
77 = ^/u2i in view of (|2.18j) and ()3.2j) . The other is between and y/u-z, in view of ()4.6j) . At 77 = 0, 
the function has a positive value. To see this, 

(4.7) -xP o (0, 1, u 2 , 0) + 80tP 2 (0, 1, u 2 , 0) = P (0, 1, u 2 , 0)[-x + t/i 3 (l, u 2 , 0)] . 
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According to Lemma 3.4, /i 2 (l, u 2 , 0) > /U 3 (l, u 2 , 0) when u 2 > u* or equivalently when a < x/t < 16. 
It follows from formula (|2.4jl and inequality (j2.11j) that Po(0, 1, u%, 0) < 0. Hence, the right hand side 
of (j4.7j) is bigger than 

P (0,l,U2,0)[-x + t/x 2 (l,u 2 ,0)] = , 

where the equality comes from (j3.2|) . Since it is a cubic polynomial in rf and since it is positive for 
large r\ > 0, the function — xPo(r] 2 , 1, m 2 , 0) + 80tP 2 (7/ 2 , 1, w 2 , 0) can have at most two zeros on the 
positive 77-axis. Hence, the above two zeros are all simple zeros. 

It now becomes straight forward to check the variational conditions (j4.2H4.3j) and the constraint 
(j4.4jl . just as we do in case (2). 

□ 

5. Other Step Like Initial Data 

In this section, we will classify all types of step like initial data (jl.7j) for equation (jl.4j) . When 
a = 0, since b ^ 0, the solution of (jl.5j) will never develop a shock. We therefore study the cases 
a > and a < 0. In the former case, it is easy to check that, when b > a, the solution of equation 
(jl.5j) will never develop a shock; accordingly, we will restrict to b < a. Similarly, in the latter case, 
we will confine ourselves to b > a. 

We will only present our proofs briefly, since they are, more or less, similar to those in Section 3. 

5.1. Type I: a > 0, a/4 < b < a. 

Theorem 5.1. (see Figure 2.) For the step like initial data with a > 0, a/4 < b < a, the 

solution of the Whitham equations HI. 6}) is given by 



Ui — a , x = /x 2 (a, m 2 , 6) t , u 3 = b 
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Figure 2. Self-Similar solution of the Whitham equations for a = 1 and b = 1/3 of 
type I. 

for /i 2 (a, b,b) < x/t < /i 2 (a, a, b), where /i 2 (ot, b, b) = —10a 2 — 40ab + 806 2 and // 2 (a, a, b) = 16a 2 + 
8ab + 66 2 . Outside this interval, the solution of M.5)) is given by 

u = a x/t < /i 2 (a, 6, 6) 

and 

u = b x/t > /U 2 (a, a, 6) . 

Proof. It suffices to show that [12(0, U2,b) is an increasing function of U2 for < w 2 < a. By (|2.28|) . 
we have 

dM(a,u 2 ,b) 10(3u 2 + b - a) 

~a = 1 K - > U 

ou 2 a — 

for b < U2 < a, where we have used a/4<6<ain the inequality. Since M(a, -u 2 , b) = at w 2 = 0, this 
implies that M(a, u 2 , 6) > for b < u 2 < a. It then follows from ()2.27|) that /i 2 (a, u 2 , 6)— /i 3 (a, u 2 , 6) > 
0. By Lemma 2.2, we conclude that 

ci/i 2 (a,M 2 ,o) 
dii 2 

for 6 < m 2 < a. 

□ 
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5.2. Type II: a > 0, -2o/3 < & < a/4. 

Theorem 3.1 is a special case of the following theorem. 

Theorem 5.2. (see Figure 1.) For the step like initial data \1.T\) with a > 0, — 2a/3 < b < a/4, the 
solution of the Whitham equations M.b\) is given by 

u 1 = a, x = fx 2 (a, u 2 , u 3 ) t , x = /j 3 (a,u 2 ,u 3 ) t 

for —15a 2 < x/t < fi 2 (a, u**, b) and by 

ui = a , x = \i 2 {a, u 2 , b) t , u 3 = b 

for fi 2 {a,u** ,b) < x/t < 16a 2 + 8ab + 6b 2 , where u** is the unique solution u 2 of n 2 (a,u 2 ,b) = 
A*3(a, u 2 ,b) in the interval b < u 2 < a. Outside the region —15a 2 < x/t < 16a 2 + 8ab + 66 2 ; the 
solution of the Burgers type equation JO} ) is given by 

u = a x/t < —15a 2 

and 

u = b x/t > 16a 2 + 8ab + 66 2 . 
Proof. The trailing edge is determined by 
(5.1) F(a,u 2 ,u 3 ) = 

when u 2 = u 3 . Here F is given by ()2.29j) . In view of the expansion ()2.3()j) . the above equation when 
u 2 = u 3 , i.e., s = 0, reduces to 

3 

2(a-u 3 ) - -(a + Au 3 ) = , 
which gives u 2 = u 3 = a/4 at the trailing edge. 
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Having located the trailing edge, we solve equation 1)5.1)) in the neighborhood of u 2 = u 3 = a/4. 
We use the expansion ()2.3()j) to calculate 

dFja, % f ) = dF(a, f , f ) = ^ 
du 2 du 3 

which implies that equation 1)5.1)) can be solved for u 3 as a decreasing function of u 2 near u 2 = u 3 = 
a/4. 

The solution of 

(5.2) /j, 2 (a,u 2 ,u 3 ) - fi 3 (a,u 2 ,u 3 ) = 

can be extended as long as a > u 2 > a/4 > u 3 > b. To see this, we need to show that 
dn 2 (a,u 2 ,u 3 ) _ dfi 3 (a,u 2 ,u 3 ) _ dfx 2 (a,u 2 ,u 3 ) dfi 3 (a,u 2 ,u 3 ) 

r% U, _ U, ^U, U 

OM3 ou 2 ou 3 

on the solution of ()5.2)) . The proof of the equalities is the same as that of (j3.9j) in Section 3. To 

prove the inequalities, in view of Lemma 2.2, it is enough to show that 

dq(a,u 2 ,u 3 ) dq(a,u 2 ,u 3 ) 
du 2 du 3 

We use formulae ()2.19j) to rewrite equation (j5.2j) as 

1 n 0/ , , ^dq(a,u 2 ,u 3 ) 1 gg(a,w 2 ,M 3 ) 
- A 2 - 2(a + u 2 + u 3 ) 5 = o A3 ~ 2 ( a + M 2 + u 3 ) ■= , 

which, together with inequalities ()2.14|) and ()2.15)) . proves that and ^ have the same sign on 
the solution of ()5.2j) . On the other hand, we calculate from ()2.22)) 

<9g(a,M 2 ,M 3 ) 

= 4(a + 3m 2 + u 3 ) > 

aw 2 

for a > u 2 > a/4 > u 3 > b > — 2a/3. 

We now extend the solution of ()5.2|) as far as possible in the region a > u 2 > a/4 > u 3 > b. There 
are two possibilities: (1) u 2 touches a before or simultaneously as u 3 reaches b and (2) u 3 touches b 
before u 2 reaches a. 
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Possibility (1) is impossible. To see this, we use ()2.17|) and ()2.19j) to calculate 
(5.3) /x 2 (a, a, u 3 ) - /i 3 (a, a, u 3 ) = 2(a - tt 3 ) 9<1 " 3 = 8(a - w 3 )(2a + 3m 3 ) , 

which, in view of b > — 2a/3, is positive for b < u 3 < a. 

Therefore, u 3 will touch b before u 2 reaches a. When this happens, we have /i 2 (a, w 2 ,o) — 
/^ 3 (a, w 2 ,o) = 0. In the same way as we prove Lemma 3.4, we can show that this equation has 
a unique solution w 2 in the interval b < w 2 < a. 

The rest of the proof is similar to that of Theorem 3.1. 

□ 









- u 
























x/t 


-30 


-20 






10 


20 
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Figure 3. Self-Similar solution of the Whitham equations for a = 1 and b = —2/3 of 
type III. 

5.3. Type III: a > 0, b = -2a/3. 

Theorem 5.3. (see Figure 3.) For the step like initial data with a > 0, b = —2a/3, the solution 
of the Whitham equations M.b\) is given by 

U\ — a , x = /i 2 (a, «2) u 3 ) t , x = fi 3 (a, u 2 , u 3 ) t 

for —15a 2 < x/t < 40a 2 /3. Outside the region, the solution of the Burgers type equation 4 1.5]) is 
given by 

■u = a x/t < —15a 2 
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and 



u = b x/t> 40a 2 /3 



Proof. It suffices to show that u 2 and U3 of ^2(0, U2, U3) — /-^(a, U2, U3) = reaches a and b = — 2a/3, 
respectively, simultaneously. To see this, we deduce from equation ()5.3|) that 



(5.4) 



// 2 (a, a, -2a/3) - /i 3 (a, a, -2a/3) = 8(a - 2a/3)[2a + 3(-2a/3)] = 



□ 
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Figure 4. Self-Similar solution of the Whitham equations for a = 1 and 6 = —1.1 of 
type IV. 

5.4. Type IV: a > 0, b < -2a/3. 

Theorem 5.4. (see Figure 4-) For the step like initial data Ji.7| ) with a > 0, b < —2a/3, the solution 
of the Whitham equations M.b]) is given by 

ui = a , x = n 2 (a, u 2 , u 3 ) t , x = 7/3 (a, u 2 , u 3 ) t 

for —15a 2 < x/t < 40a 2 /3. Outside the region, the solution of the Burgers type equation M.5)) is 
given by 

u = a x/t < —15a 2 

and 



u 



s /^£ 40a 2 /3 < x/t < 306 2 
b x/t > 306 2 
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Proof. By the calculation (|5.4|) . when u 2 of 712(0, u 2 , u 3 ) — /i 3 (a, xx 2 , M3) = touches a, the correspond- 
ing U3 reaches — 2a/3, which is above 6. Hence, equations 

x = 7x2(0, U2, U3) t , x = n 3 (a,u 2 ,u 3 ) t 

can be inverted to give 1*2 and XX3 as functions of x/t in the region /x 2 (a, a/4, a/4) < x/t < 7x2(0, a, — 2a/3). 
To the right of this region, the Burgers type equation (J1.5)) has a rarefaction wave solution. 

□ 



x/t 
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Figure 5. Self-Similar solution of the Whitham equations for a = — 1 and 6 = of 
type V. 

5.5. Type V: a < 0, 6 < -a/4. 

Theorem 5.5. (see Figure 5.) For the step like initial data \1.7\) with a < 0, a < b < — a/ A, the 
solution of the Whitham equations M.b]) is given by 

U\ = b , x = 7x2(0, 1*2; a) t , % = a 

for 7x2(0, 0, a) < x/t < 7x2(0, a, a), where 7x2(6, 6, a) = 6a 2 + 8a6 + 166 2 ana 1 7x2(6, a, a) = 80a 2 — 40a6 — 
106 2 . Outside this interval, the solution of is given by 

u = a x/t < 7x2(6, 6, a) 
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and 

u = b x/t > fi 2 (b,a,a) . 

Proof. It suffices to show that « 2 (a, w 2 ,6) is a decreasing function of u 2 for a < u 2 < b. By (|2.19|) . 
we have 

dfi 2 (b,u 2 ,a) ld\ 2 dq 1 <9 2 g 

The second term is negative because of (|2.14j) and = 12 > 0. The first term is also negative. Its 
first factor is positive in view of Ijl.lOj) . The second factor 

dq 



du 2 

for a < u 2 < b because of b < —a/4. 



4(6 + 3u 2 + a) < 



□ 
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Figure 6. Self-Similar solution of the Whitham equations for a = — 1 and b = 1.2 of 
type VI. 

5.6. Type VI: a < 0, -a/4 < b < -2a. 

Theorem 5.6. (see Figure 6.) For the step like initial data \1. 7\ ) with a < 0, —a/4 < b < —2a, the 
solution of the Whitham equations M.b]) is given by 

(5.5) x = /ii(wi,u 2 ,a) t , x = j^ 2 (ui,u 2 ,a) t , u^ = a 
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for ha 2 < x/t < fi 2 (b, u***, a) and by 

(5.6) ui = b , x = /i 2 (b,u 2 ,a) t , u 3 = a 

for fi 2 (b, u***, a) < x/t < 80a 2 — AOab — 10b 2 , where u*** is the unique solution u 2 of fii(b,U2,a) = 
H2{b,u 2 ,a) in the interval a <u 2 <b. Outside the region 5a 2 < x/t < 80a 2 — A0ab — 10b 2 , the solution 
of the Burgers type equation fPJ) is given by 

u = a x/t < 5a 2 

and 

u = b x/t> 80a 2 - AOab - 10b 2 . 

Proof. We first locate the "leading" edge, i.e., the solution of equation (|5.5|) at u\ = u 2 . Eliminating 
x/t from the first two equations of (j5.5J) yields 

(5.7) /ii(ui,u 2 ,a) - /x 2 (mi,M2,o) = . 

Since it degenerates at Ui = u 2 , we replace (|5.7|) by 

/p- q\ nl \ Vi{ui,u 2 ,a) - n 2 {u x ,u 2 ,a) 

(5.8) G{u x , u 2 , a) := — — = . 

(mi - u 2 )K{s) 

Using formulae ()2.19j) for fix and fi 2 and formulae (|2.12j) for Ai and X 2 , we write 

- E[E -l- a)K] ^ + a lt} E ~ i 1 - ^Z K} ■ 

In view of (12. 7j) and (|2.8|l . equation ()5.8j) reduces to 

<9g(ui, w 2 , a) dq(ui, u 2 , a) _ n 



dui du 
at the "leading" edge u\ — u 2 . This gives 



a 

Ui = u 2 = -- 
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Having located the "leading" edge, we solve equation (|5.8|) near u\ = u 2 = —a/4. We calculate 

dG(-a/4,-a/4, a) _ dG(-a/4,-a/4,a) _ ^ 
du\ du 2 

These show that equation (|5.8jl gives U\ as a decreasing function of u 2 

(5.9) mi = B(u 2 ) 

in a neighborhood of U\ = u 2 = —a/4. 

We now extend the solution ()5.9|) of equation ()5.7|) as far as possible in the region a < u 2 < —a/4 < 
Ui < b. We use formula (j2.19|) to calculate 

dfix ldXx dq 1 <9 2 g 

dux 2duidui 2 9«j 

<9/i 2 1<9A 2 <9g 1 <9 2 g 
= o^— a h o A 2 - 2(«i + u 2 + a) — T . 

In view of ffTTTUjl . (|2~T3j) and (EH, we have 

^ > if |«->0, 
oui au\ 

— - < if — < . 

ou 2 ou 2 

We claim that 

(5.10) p->o, |^<0 

on the solution of (|5.7|) in the region a <u 2 < —a/4 <U\ <b. To see this, we use formula (|2.19|) to 
rewrite equation ()5.7|) as 

1 dq 1 (9a 

-[Ai - 2(«i + u 2 + a)]— i- = -[A 2 - 2( Ul + m 2 + a) 1 ' 



2 du\ 2 <9ii 2 

This, together with 

* S « 2 ( Ml -„ 2 )^ = 8 ( Bl -„ 2 )>0 



<9wi du 2 du\du 2 
for «i > u 2 , and inequalities ()2.13|) and ([2. 14)1 . proves ()5.10|) . 
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Hence, the solution (|5.9|) can be extended as long as a < u 2 < —a/4 < u\ < b. 
There are two possibilities: (1) U\ touches b before u 2 reaches a and (2) u 2 touches a before or 
simultaneously as u\ reaches a. 

Possibility (2) is impossible. To see this, we use (I2.16J) . ()2.19|) and ()2.22|) to calculate 

(5.11) ^i(wi) a, a) — /-^(wi, a, a) = 40(«i — a)(ui + 2a) , 

which is negative for —a/4 < u\ < b < —2a. 

Therefore, u\ will touch b before u 2 reaches a. When this happens, we have 

(5.12) fii(b, u 2 ,a) — // 2 (&, u 2 ,a) = . 

Lemma 5.7. Equation has a simple zero, counting multiplicities, in the interval a < u 2 < b. 

Denoting this zero by u*** , then fi 1 (b,u 2 ,a) — n 2 (b,u 2 ,a) is positive for u 2 > u*** and negative for 
u 2 < u*** . 

Proof. We write 

/ \ /7 \ n \ 2(b — u 2 )K dq dq^^ . . dq . 

(5.13) ^(6, u 2 , a) - ^ 2 (b, u 2 , a) = E[E _ {1 _ s)K] {(^ + s-^)E - (1 - s)—K} . 

Denote the parenthesis of (|5.13j) by N(b, u 2 , a). Since E — (1 — s)K > for a < u 2 < b, the left hand 
side has a zero iff N(b, u 2 , a) on the right has one. 
We now calculate 

dN(b,u 2 ,a) _ 30E(s) a-b 
du 2 b — a 3 

Since N(b, u 2 , a) is zero at u 2 = a and positive at u 2 = 6, we conclude from the above derivative that 
N(b, u 2 , a) has a simple zero in a < u 2 < b. □ 

We now continue to prove Theorem 5.6. Having solved equation (J5.7|) for u\ as a decreasing 
function of u 2 for u*** < u 2 < —a/4, we can then use the last two equations of ()5.5j) to determine U\ 
and u 2 as functions of x/t in the interval fi 2 (—a/4, —a/4, a) < x/t < fi 2 (b, u***, a). 
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We finally turn to equations (j5.6jl . We want to solve the second equation of (|5.6|h x/t = ^{b, u 2 , a), 
for u 2 < u*** . It is enough to show that /i 2 (6, u 2 , a) is a decreasing function of u 2 for u 2 < u***. 

According to Lemma 5.7, /ii(6, w 2 ,a) — /i 2 (6, w 2 ,a) < for u 2 < «***. Using formula (j2.19|) for //i 
and /z 2 , we have 

I [Al _ 2(6 + u 2 + a)]|^- < ±[A 2 - 2(6 + u 2 + a)]p- . 

2 OUi 2 OU 2 

This, together with 

9q dq 2(b-u 2 )^ = 8(b-u 2 )>0 



du\ du 2 du\du 2 
for u\ > u 2 , and inequalities ()2.13|) and ()2.14|) . proves 

dq(b, u 2 , a) 



du 2 



< 



for u 2 < u*** . Hence, 



du 2 ld\ 2 dq l r , „ ^d 2 q 

£r-2i^8k + 2 lx '- 2ib+U2+a)] ai <0 



□ 




5.7. Type VII: a < 0, b = 



-2a. 
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Theorem 5.8. (see Figure 7.) For the step like initial data \1. 7\ ) with a < 0, b = —2a, the solution 
of the Whitham equations M.n\) is given by 

x = fi\(ui, U2, a) t , x = fiz(ui, U2, a) t ,Uz = a 

for 5a 2 < x/t < 120a 2 . Outside the region, the solution of the Burgers type equation M.b^l is given 
by 

u = a x/t < 5a 2 

and 

u = b x/t> 120a 2 . 

Proof. It suffices to show that u\ and u<i of /ii(wi,M2,a) — fi2(ui,U2,a) = reaches b = —2a and a, 
respectively, simultaneously. To see this, we deduce from equation (|5.1H1 that 

(5.14) fJ>i{ui, a, a) — /^(/Ui, a, a) = 8(ui — a){u\ + 2a) 

is negative for u\ < b and vanish when u± = b = —2a. 

□ 

5.8. Type VIII: a < 0, b > -2a. 

Theorem 5.9. (see Figure 8.) For the step like initial data \1. 7| ) with a < 0, b > —2a, the solution 
of the Whitham equations M.b]) is given by 

x = fii(ui, U2, a) t , x = /^(wi, «2, a) t , u^ = a 

for 5a 2 < x/t < 120a 2 . Outside the region, the solution of the Burgers type equation il.ty) is given 
by 

■u = a x/t < 5a 2 
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and 

{yf^ t 120a 2 < x/t < 306 2 
■ 
b x/t > 306 2 

Proof. By the calculation (|5.14|) . when u 2 of fi 2 (ui,u 2 , a) — fi 3 (ui,U2,a) = touches a, the corre- 
sponding u 3 reaches —2a, which is below 6. Hence, equations 

x = fi 2 (a,u 2 ,u 3 ) t , x = fi 3 (a,u 2 ,u 3 ) t 

can be inverted to give u 2 and u 3 as functions of x/t in the region fi 2 (— a/4, — a/4, a) < x/t < 
fi 2 (—2a, —2a, a). To the right of this region, the Burgers type equation (jl.5j) has a rarefaction wave 
solution. 

□ 




Figure 8. Self-Similar solution of the Whitham equations for a = —1 and b = 2.5 of 
type VIII. 
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